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A comprehensive and detailed analysis of hadronic abundances measured in Au-Au collisions at 
RH1C at \/s NN = 130 and 200 GeV is presented. The rapidity densities measured in the central 
rapidity region have been fitted to the statistical hadronization model and the chemical freeze-out 
parameters determined as a function of centrality, using data from experiments BRAHMS, PHENIX 
and STAR. The chemical freeze-out temperature turns out to be independent of centrality to a few 
percent accuracy, whereas the strangeness under-saturation parameter 7s decreases from almost 
unity in central collisions to a significantly lower value in peripheral collisions. Our results are in 
essential agreement with previous analyses, with the exception that fit quality at \/s NN = 200 GeV 
is not as good as previously found. From the comparison of the two different energies, we conclude 
that the difference in fit quality, as described by \ 2 values, is owing to the improved resolution 
j /' of measurements which has probably exceeded the intrinsic accuracy of the simplified theoretical 

formula used in the fits. 

O 
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I. INTRODUCTION 

The idea of a statistical model to account for multiple hadron production processes in high energy collisions dates 
■ back to a work by Fermi [l| . This model has been successful in reproducing the production rates of measured 
hadronic species in collisions of both elementary particles @, [H, H[ and heavy-ions @, [y, 0, H, @ • These models have 
been extensively and succesfully applied to the phenomenon of multifragmentation in nuclear collisions (loj . 

Many evidences have been collected that a new form of matter, the Quark-Gluon Plasma (QGP), where effective 
degrees of freedom are quarks and gluons, has been created in the collisions of heavy ions at relativistic energies at 
[> ■ top SPS and RHIC energy. While the success of the statistical-thermal model in describing particle multiplicities 
indicates that local thermodynamical equilibrium has been achieved and seems to confirm the general hypothesis of 
• QGP formation, some aspects are still to be understood, like e.g. whether strangeness is fully equilibrated and the 
relation of this successful description with that in elementary collisions. In this respect, an analysis of RHIC data, at 
a center-of-mass energies 130 and 200 GeV per collid ing pair of nucleon, can be illuminating. Indeed, similar analyses 
have been carried out in the past few years [ill, EH fl3L Il4 EH E3|- Yet, recently, a bunch of new experimental 
data has been published which makes it a worthwhile step to provide an independent analysis including this newly 
available data. Furthermore, we have studied hadron pro duction at lower beam energies in heavy-ion collisions at 
AGS and SPS with the statistical hadronization model [l8l . Ell and with this paper we complete our previous works. 

There is also another issue which motivates our analysis. So far, most analyses have been using as an input to 
the fit (N — 1) particle ratios formed a posteriori from N measured particle multiplicities without including any 
normalizing yield. This procedure was based on the tacit assumption that fitting either ratios or multiplicities lead 
to equivalent results. This is in general not true and the outcome of a statistical analysis relying on only ratios of 
hadron multiplicities may be seriously biased po| , depending on the input set of ratios. The reason of this problem 
is that, in principle, one can form N(N — l)/2 different combinations 1 of particle ratios from N different measured 
multiplicities and choosing a subset of them implies an information loss. Moreover, the different ratios are obviously 
correlated if a particle appears more than once in the ratios and those correlations must be taken into account in 
the x 2 minimization, thus complicating the fit. In the worst case, the central values of the fitted statistical model 
parameters may deviate several standard errors from the central values of parameters determined from a fit to particle 
multiplicities whereas the actual magnitude of the error is not possible to know without explicit comparison. We stress 
that this problem arises when using ratios calculated a posteriori from a set of primordially measured multiplicities, 
while ratios directly measured by the experiments because of beneficial systematic error cancellation (e.g. n~ /n + or 
p/p) are perfectly safe. 

While at lower beam energies integrated multiplicities in full phase space are a more suitable input for the statistical 
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model at RHIC energies of v^nn > 1^0 GeV, rapidity distributions are wide enough to allow the extraction of 
the thermodynamical properties of the average fireball produced at mid-rapidity with rapidity densities themselves. 
In fact, the standard width of charged particle rapidity distribution at -\/%n = 200 GeV is 2.1 [2l[, sufficiently 
larger than standard widths of single-fireball rapidity distribution (at most 0.8 for pions at the kinetic freeze-out 
temperature of 125 MeV). On the other hand, at top SPS energy v^nn = 1^.2 GeV the measured rapidity width is 
1.3 [22l | . which is consistently smaller and closer to the single-fireball width. Therefore, using mid-rapidity densities 
at this and lower energies artificially enhances heavier particles with respect to lighter ones as they have a narrower 
rapidity distribution. 



II. THE DATA ANALYSIS 



We have analyzed the rapidity densities in Au-Au collisions at 130 and 200 GeV per participating nucleon measured 
by BRAHMS, PHENIX and STAR collaborations at RHIC employing a version of statistical hadronization model 
described in detail in [l8l . [l9j . The use of the grand-canonical formalism is appropriate here in that particle multi- 
plicities are large. The formula for the ith primary hadron (including both stable hadrons and resonances) rapidity 
density reads: 



drij 
dy 



dV (2Ji + 1) 
dy (2tt) 3 



d 3 p 



7s 



-n s i e V , P 2 +mf/T-/x-q i /T ± 1 



(1) 



where T is the temperature, = (Qi, Bi, Si) is a vector having as components the electric charge, baryon number 
and strangeness of the hadron and fj, — (fJ,Q, fJ-B, A*s) is a vector of the corresponding chemical potentials; 75 is the 
strangeness under-saturation factor and n s i is the number of valence strange quarks in the ith hadron; the upper sign 
applies to fermions, the lower to bosons. The absolute normalization dV/dy in {1} is a product of the rapidity density 
of clusters at mid-rapidity p(0) times the volume of the average fireball at mid-rapidity [l9| . For the above formula to 
make sense, the parameters T, and strangeness under-saturation factor 75 should be constant over a rapidity range 
encompassing the single fireball rapidity width jl9| |. In this work, the chemical potential us is determined enforcing 
vanishing strangeness density and /iq by requiring the final ratio of charge to baryon number to equalize the initial 
one Z/A, i.e. by assuming that there is no major dependence of the strangeness and electric density on rapidity. The 
other 4 parameters (T, [is, 7s, dV/dy) are determined by minimizing the x 2 : 



x 2 = 



E 



(dNf/dy - dNl/dy) 2 



{dNf/dy - n\ dVjdyf 



(2) 



in which dNf/dy and dNf/dy are the experimental and theoretical rapidity densities, n\ is the particle density 
evaluated within statistical hadronization model and Ci is the experimental error of the rapidity density of a particle 
species i. Unless otherwise stated, all experimental errors quoted in this paper are a quadratic sum of statistical 
and systematic errors. Before going deeper into the data analysis and exploration of results, we need to discuss a 
preliminary treatment of the experimental data which was necessary to make a combined analysis of all experiments. 



A. Centrality interpolations 

Particle rapidity densities are measured in some selected centrality windows which are different for different exper- 
iments. As a consequence, measurements relevant to e.g. the most central collisions from different experiments must 
be renormalized. Moreover, the chosen centrality windows can be different for different particle species even within 
the same experiment. Therefore, in order to make a correct analysis of the full data set, one needs to find a proper 
method to estimate rapidity densities of different hadronic species in the same centrality window, i.e. a proper inter- 
polation method. This can be done in many ways. For example, taking advantage of a possible linear and parabolic 
scaling with number of participants (Np) or number of binary collisions (Nbi n ) as well as with the negative hadron 
pseudo-rapidity density (dN h - /drj). In many cases such a simple scaling behaviour can be found. For example, STAR 
collaboration has found out [23|] that L and A rapidity densities scale well with the dN h -/dr? in Au-Au collisions at 
130A GeV. However, none of the above-mentioned scaling variables with simple functional forms is able to describe 
the centrality dependence of all different hadron species; instead, more complex dynamical combination (e.g. [24j, |25j) 
of "hard" and "soft" physics processes must be considered. 

Typically the interpolation correction is not very large. For instance, we might want to estimate the rapidity density 
of a hadron species in the [0-6%] most central collisions while the experimental value is given for the [0-5%] and [5-10%] 
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most central collisions. It is also important to note that centrality fractions are independent of the observable used to 
define them provided that the observable varies monotonically; this ensures the one-to-one correspondence between 
different observables. We have chosen an interpolation method we deem is more robust and model-independent than 
any simple scaling with Np, Nbi n or dN h - /drj and implemented it consistently for all 2 rapidity densities that need to 
be interpolated. We write the rapidity density of a hadron species i as a fc max th order polynomial of the centrality as 
follows: 



d-'.V, '' 



E a ^ fe ; c=i- x (3) 



dydc 

tf fc=0 

in which a\ denote free parameters and x G [0, 1] is the fraction of the differential cross section as a function of the 
variable defining centrality itself (0=0% most central and 1=100% most central collisions). The rapidity density of 
a hadron species in a certain centrality window [y m - m , y ma x] is obtained by integrating Eq. ([3]) over [y m i n , 2/max] ■ We 
thereby calculate the rapidity densities in the centrality windows where they have been measured and fit the fc max +l 
free parameters to reproduce the measured rapidity densities. Once the parameters are fitted, rapidity densities of 
any hadron species can be estimated in any centrality window [y' m i n , y' max \ by simply integrating with respect to c 
over the desired region of centrality keeping the fitted parameters fixed. The maximal order of the polynomial (i.e. 
maximal number of free parameters minus one) for a certain particle species i is the number of centrality bins in 
which it is measured. We have always chosen the maximal order k max =N W i n dows-l for these interpolations unless 3 
this would lead to non-monotonic behaviour of Eq. ([3]) within the centrality range we are interested in, in which case 
we have chosen the maximal order such that the polynomial is monotonically increasing. The original experimental 
errors are properly propagated so that the errors of our interpolated rapidity densities include the original errors as 
well as the additional uncertainty arising from the interpolation. 

An example of our interpolations is shown in fig. [TJ In the left panel, L (open round symbols) and A (open square 
symbols) rapidity densities measured by STAR collaboration at 130A GeV are shown while the filled symbols denote 
our estimates for the corresponding rapidity densities in 7 other centrality windows, namely the ones in which kaons 
and nucleons are measured. In the right panel of fig. [TJ similar plot is shown for S~ and S + at 200A GeV. 

The experimental rapidity densities as well as our interpolated values are shown in tables III [TT1 and ITTT1 The numbers 
in plain text are our estimates while numbers written in bold case are the experimental values our interpolations are 
based on. 



B. Au-Au collisions at 130 A GeV 

STAR collaboration has measured p [26[, p (26j, K + [27j and K~ [27| rapidity densities around mid-rapidity in 
Au-Au collisions at 130A GeV in 8 centrality bins, see table [TJ Proton and anti-proton rapidity densities include 
weak feeding from ( mult i) strange hyperons. We have chosen to perform our statistical model analysis as a function 
of centrality at this beam energy in the centrality windows nucleons and kaons are measured. 

On the other hand, A and A rapidity densities are measured in 5 different centrality bins [23| and the rapidity 
densities also include weak feeding from multi-strange hyperons. By fitting different functions to the transverse mass 
spectra, STAR collaboration has obtained two slightly different values for both A and A rapidity densities, thus we 
have taken a weighted average of those. Since A and A's centrality bins differ somewhat from the centrality bins in 
which nucleons and kaons are measured, we have estimated the hyperon yields in the latter ones as described in the 
previous section. 

For the fit to be reliable, it is necessary to include the most abundantly produced particles, charged pions, in the 
analysis. However, the integrated STAR tt^ rapidity densities (corrected for weak decay feeding) are publicly available 
only in the 0-5% most central bin at 130A GeV |28j | and so polynomial interpolation for peripheral collisions can not 
be implemented. Without better knowledge, we have assumed that the ratio 

^/^- = ^-/^- = 239.0/296.6 (4) 
dy I dr\ dy I dij 

retains its value in every centrality bin. Since most of the charged hadrons emitted in a heavy-ion collision are pions, 
we deem that Eq. ([4| yields a reasonable estimate of the pion rapidity densities. As far as the error on the estimated 



2 there are few exceptions to this rule which will be explicitly discussed 
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FIG. 1: Left panel (a): Experimental L (open round symbols) and A (open square symbols) rapidity densities as a function 
of centrality (c m in and c max being the borders of the centrality window under consideration) in Au-Au collisions at 130A 
GeV measured by STAR collaboration [23j]. Full symbols denote our interpolated rapidity densities in centrality windows the 
statistical analysis is carried out. 

Right panel (b): Measured E~ (open round symbols) and H + (open square symbols) rapidity densities in Au-Au collisions at 
200A GeV [17(. The full symbols denote our estimates for the same rapidity densities in the centrality windows the statistical 
analysis is carried out. The round symbols are shifted 0.02 leftward for clarity. 

c m in and c m ax denote the limits of the centrality bin corresponding to each data point; in other words, the data points have 
been stick in the centre of the bin. 



rapidity densities is concerned, we have have added the additional systematic error arising from the extrapolation to 
different centrality bins in quadrature with the relative error of 4.5% that is quoted in the 0-5% most central collisions. 
Based on the published [23|, [2?| systematic errors of dN h -/d?7 in the various centrality bins, we estimate an overall 
error of 10.3% in the pion rapidity density. In order to check the stability of our analysis, we have repeated the fits 
by assuming 5% and 15% error in the pion rapidity densities in each centrality bin. The ensuing fitted parameters 
showed little difference and their central values turned out to be well within the error bar of the main fit. 

The rapidity densities of hyperons S~ and S + have been measured in 3 centrality windows (0-10%, 10-25% and 

25- 75% most central collisions) while 51 + Ct is measured in the 0-10% most central collisions [29j only. Similarly to 
L and A , two slightly different rapidity densities are quoted for both S's and we have taken the weighted average 
as our input for the analysis. Finally, ^-meson is measured in 3 different centrality windows (0-11%, 11-26% and 

26- 85% most central collisions) [30]. Similarly to A's, we have estimated the and </> rapidity densities in the 8 
STAR reference centrality bins, but since data is available in 3 centrality windows only, a 2nd order interpolation 
polynomial was used in Eq. ([3]). 

Because of the vast width of the peripheral centrality window (w 30 — 80% most central collisions), our method 
to estimate S and (f> rapidity densities fails in the most peripheral bin (58-85%), and in general the relative errors 
increase with decreasing centrality. Particularly, the extrapolation of Q + Cl rapidity density from central to peripheral 
collisions based on a single centrality is meaningless. Thus, we have removed ft from the STAR particle set to estimate 
the freeze-out parameters in Au-Au collisions at 13QA GeV in the 7 most central bins. As a check, the analysis has 
been repeated in the two most central bins by including + fl rapidity density, which we have assumed to scale with 
the negative hadron pseudo-rapidity density, at least in this short range. The fit outcome turns out to be essentially 
unaffected by this inclusion and thus all quoted results in this work at 130A GeV refer to fits without fl's. 

Also PHENIX collaboration has measured 7r + , 7r~, K + , K~ , p and p rapidity densities in a pseudo-rapidity window 
of 1 77 1 < 0.35 around mid-rapidity in Au-Au collisions at 130 A GeV [3l|,[33]. The data is divided in 5 centrality bins 
(see table U) which differ from STAR's ones and thus direct comparison is not possible. Also, A and A rapidity 
densities have been measured in the most central bin at 130A GeV (33| . No weak decay corrections were applied to 
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any of the hadron species. We have repeated the fits with the PHENIX data and found out that the data set is rich 
enough to fix all the statistical model free parameters only in the most central bin, in which hyperons are included in 
the data sample (see fig. [2] and table IIV| ) . In the other bins, the set n^, if , p and p does not allow to reliably fit all 
the 4 free parameters because of the relatively short lever arm in mass and the low baryon chemical potential which 
makes pions multiplicities too close. 

Finally, we have made a combined fit to PHENIX and STAR data. First, the PHENIX rapidity densities of n ± , 
K^, p and p have been estimated in the STAR centrality bins (see again table [TJ) according to the aforementioned 
interpolation procedure. The obtained and nucleon rapidity densities agree very well with the corresponding 
experimental STAR values in the most central bin while the relative discrepancy increases in the more peripheral ones, 
yet within the error bars. It should also be pointed out that PHENIX rapidity densities of pions are larger than the 
corresponding STAR values because in the former case no weak decay corrections was applied. The possible different 
overall normalization between the 2 experiments was taken into account by introducing one more free parameter fp 
multiplying all PHENIX rapidity densities; otherwise stated, the common scaling factor dV/dy becomes fp x dV/dy 
for PHENIX data. Obviously, one expects fp » 1 in each centrality bin, if the experiments are to be in essential 
agreement. 

The statistical model best fit parameters determined from the combined STAR+PHENIX fit are shown in fig. [2] 
and table IIVI along with those determined by a fit to STAR data alone. It can be seen that they are in very good 
agreement with each other and that the cross-normalization parameter fp varies between 0.8 to 0.92 throughout the 
examined centrality range. 



Au-Au collisions at 200A GeV 



The analysis has been carried out also at v^nn = ^00 GeV. At this energy, PHENIX collaboration has measured 34] 
the rapidity densities of the same set of hadrons as at \/s NN = 130 GeV, i.e. ir + , ir~ , K + , K~, p and p over 11 
centrality bins. Pion rapidity densities do not include weak decay products while proton and anti-proton rapidity 
densities are corrected from A and A feeding. Like at the lower beam energy, the set of different hadron species is not 
large enough for us to fit the statistical model parameters reliably. 

STAR collaboration has measured the same hadron species (n^, K ± , p, p, A, A, (f>, S T and fi + fi E3, HE) H(| ) as 
at V*nn = 130 GeV. The centrality windows in which STAR and PHENIX have measured pions, kaons and protons 
are mostly the same. In order to minimize the amount of data manipulation, we have chosen to keep the PHENIX 
data as it is and to estimate STAR rapidity densities, whenever necessary, in PHENIX centrality bins. The hyperons 
and (f> meson have been measured in many more centrality windows than at the lower beam energy. Still, most of the 
windows are wider than the ones for pions and so we have interpolated hyperon rapidity densities in the narrower 
pions centrality windows according to Eq. (J3]). All rapidity densities measured by STAR collaboration are cleaned 
from weak decay products except p and p which include feeding from hyperons. 

Finally, BRAHMS collaboration has measured the same hadron rapidity densities [13] as PHENIX collaboration in 
4 different centrality windows. The pion rapidity densities do not include any weak decay products while only L and 
A decay products are subtracted from nucleon rapidity densities. We have estimated the BRAHMS rapidity densities 
in the same centrality bins defined by PHENIX collaboration. 

Similarly to the lower beam energy, we have determined the chemical freeze-out parameters by performing a fit to 
STAR data alone and then a combined fit to STAR, PHENIX and BRAHMS rapidity densities. In the combined fit, 
free parameters fp and fp multiplying theoretical rapidity densities of PHENIX and BRAHMS respectively have been 
introduced in order to take into account possible discrepancy in overall normalization among different experiments. 
As it seems that there is a significant discrepancy in the A/p among the three experiments, we have decided to exclude 
proton and anti-proton rapidity densities measured by PHENIX and BRAHMS collaborations in the analysis (see 
detailed discussion in Section iHi]) . The resulting statistical model best fit parameters are shown in table IVl and fig. [3] 



D. Further notes 



We have left out from our analysis some additional rapidity densities of hadron species that are measured at RHIC. 
The PHENIX measurement [HI of <j> meson at 200A GeV is left out from our analysis due to the very large discrepancy 
with the corresponding STAR values. We have compared the statistical model predictions for <fi meson production 
with the PHENIX measurement though, and found out a severe disagreement between the statistical model prediction 
and the PHENIX measurement (of the same order as between STAR and PHENIX measurements) . 

Also, STAR measurements [3!l[4(| of strange resonances if (892), E(1385) and A(1520) are left out from our analysis 
due to their very short lifetime which makes their decay products rescatter after chemical freeze-out, a known issue 
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(see e.g. [181 ]) in statistical model analysis in heavy-ion collisions. On top of above measurements, we have omitted 
the STAR Kg [27| measurement in our analysis. Within the statistical hadronization model, K s multiplicity is always 
between the K + and K~ yields while the STAR measurement suggests much smaller Kg multiplicity compared with 
both K + and K~ . Thus, in this case statistical hadronization model would not be able to reproduce all K rapidity 
densities on a satisfactory level and so we have decided to rely on the K yields only. In order to take into account 
the additional uncertainty on parameters implied in fits with x 2 /dof > 1, parameter errors have been rescaled by 
\/ x 2 /dof if this is larger than 1, according to Particle Data Group procedure [ill ]. 

III. DISCUSSION 

Looking at the figs. [2] and [3[ the most striking feature of statistical hadronization model fits is that temperature and 
baryon-chemical potential do not show much dependence on centrality. Particularly, temperature is constant at few 
percent level. The strangeness phase-space under-saturation parameter 75 seems to be somewhat smaller than unity 
in peripheral collisions but reaches unity in semi-central collisions and then saturates. All of this is in agreement with 
previous findings [3, [13] • Indeed, the increasing trend of 75 as a function of centrality is more evident at 130 GeV 
than at 200 GeV; furthermore, at 130 GeV 75 apparently exceeds 1 in the most central collisions. However, given 
the large error bar, this parameter is still consistent with its natural saturation value, i.e. 1. With the present level 
of accuracy, we believe that no claim can be made about different values of 75 at the two energies. 

Comparing the statistical hadronization model parameters among the two different beam energies with the same 
Np, we see very little differences. The resulting chemical freeze-out temperatures, 7s factors and scaling volumes are 
very similar and we can see mild beam energy dependence in the baryon chemical potential only. Thus, it seems that 
at RHIC energies we have reached a saturation limit in which hadrons decouple from the strongly interacting system 
at mid-rapidity in almost the same thermodynamical state. We then easily predict, in agreement with others, that 
Pb-Pb collisions at the LHC will find T w 170 MeV and 75 w 1. 

The statistical hadronization model describes the STAR data very well at 130A GeV in every centrality bin. The 
X 2 /dof is less than 1 (see table IIVI and fig. [2]) both with STAR data alone as well as with the combined STAR + 
PHENIX data. The resulting cross normalization factors fp in the combined fit are around 0.9 in the central and semi 
central collisions while in the peripheral systems we find rather low factors of the order of 0.8. The same tendency 
is already visible when comparing the experimental STAR rapidity densities and our interpolated PHENIX rapidity 
densities at 130A GeV. Up to what extent this is a manifestation of a true difference in absolute normalization among 
these two experiments in peripheral collisions or a fit artefact is not possible to decide based on the published data 
available because there is no overlap in the particle sample. At this energy, it seems that statistical hadronization 
model tends to overestimate the proton and anti-proton rapidity densities while other particle species are very well 
described and no systematic discrepancy is seen between data and model. 

Conversely, at \/s NN = 

200 GeV, the x 2 /dof values (see tableland fig-EJ), are larger than one. This is simply due 
to the better accuracy of measurements at the larger beam energy while the average relative deviations of calculated 
rapidity densities from the experimental ones are of the same order at both energies (in fact, they are 9.7% and 11.5% 
in central heavy- ion collisions at 130^4 GeV and 200A GeV respectively). The residuals (defined as the ratio between 
the difference model-data and the experimental error) and relative deviations of measured and calculated rapidity 
densities in central Au-Au collisions at 130 and 200/1 GeV are shown in tables PVTl and IVUl where one can see that the 
model is able to reproduce the data at the same level of accuracy at both beam energies. 

For this reason, it should be stressed that the \ 2 test should be used very carefully in order to avoid naive and hasty 
judgements about the validity of the model. Indeed, what these values at different energies tell us is that the simple 
formula (fTJ) is valid up to some level of resolution and fails when the accuracy of measurements exceeds it, i.e. at 200A 
GeV. This is really no surprise because formula fl} relies on several side-assumptions and approximations that are not 
expected to be exactly fulfilled. In other words, the theoretical model expressed by (fT]) is to be taken as a zero-order 
approximation and not as a precise representation of the real process. When the resolving power of experiments is 
good enough, higher order corrections become necessary, although they are very difficult to estimate and implement. 
For instance, an assumption which may not be exactly true is the vanishing of strangeness density at mid-rapidity, 
which has been used in our fits; clearly, treating /is as a further free parameter could reduce the \ 2 - Another important 
approximation is concerned with the hadron-resonance gas model, where both hadrons and resonances are handled 
as free particles with distributed mass and the contribution of non-resonating interactions among stable hadrons is 
neglected; it is clear that this approximation will fail at some very good resolution. Finally, it should be reminded 
that the sharp separation between chemical and kinetic freeze-out is also an idealization. Even though hadronic 
rescattering does not play a major role in determining particle abundances (one good evidence is the success of the 
statistical model itself), we know that it is there; thus, different inelastic reactions may cease at different stages of 
the post-hadronization expansion and this involves deviations from the simple scheme of elastic-inelastic separation. 
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FIG. 2: Chemical freeze-out temperature (top left panel (a)), strangeness under-saturation parameter 7s (top right panel 
(b)), baryon chemical potential (bottom left panel (c)) and the best fit \ 2 P er degrees of freedom (bottom right panel (d)) as 
a function of collision centrality in Au-Au collisions around mid-rapidity at 130^4 GeV. Open symbols represent fits to STAR 
data while square full symbols represent fits to combined STAR+PHENIX data and the full round symbol a fit to PHENIX 
data alone. The full square symbols are shifted 10 units rightward and the full round symbol 10 units leftward for clarity. 



Nevertheless, the x 2 fit is a useful tool to determine the best parameters of the zero-order theory but should be 
used with care as an absolute measure of the fit quality. For example, the relative errors of hadron multiplicities 
in e + e~ experiments are typically few percent only, which leads to relatively large x 2 /do/ values 01 a t least when 
compared with the x 2 /dof values in heavy- ion collisions in which the relative errors of multiplicities are typically larger 
and thus a blind comparison of the \ 2 /dof values arising from the fits to elementary collisions and heavy-ion collisions 
could be highly misleading. However, as has been mentioned, if fits have a low quality, the estimated parameter errors 
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FIG. 3: Chemical freeze-out temperature (top left panel (a)), strangeness under-saturation parameter 7s (top right panel 
(b)), baryon chemical potential (bottom left panel (c)) and the best fit \ 2 P er degrees of freedom (bottom right panel (d)) as 
a function of collision centrality in Au-Au collisions around mid-rapidity at 200A GeV. Full symbols represent fits to combined 
STAR+PHENIX+BRAHMS data while the open symbols represent fits to STAR data only. The filled symbols are shifted 10 
units rightward for clarity. 



could be unrealistically small and this is why we rescaled errors by \/x 2 /dof, according to the procedure adopted by 
Particle Data Group in such cases [ill ]. 

Fits to STAR data at 200^4 GeV do not reveal any particular discrepancy between the experiment and the model, 
all particle species are roughly equally well described (~ la — 2a deviation) with the exception of pions which are 
always very well reproduced (see fig. 0J. We have repeated the fits by systematically removing different particle 
species from the fits and found out that similarly to the lower beam energy, statistical model seems to have problems 
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here as the ratio between the difference data-model and the experimental error) distribution. The data is from STAR, the 
model values refer to the combined fit to BRAHMS, STAR, PHENIX experiments. 



in reproducing the proton and anti-proton rapidity densities especially together with the A and A rapidity densities. 
Removal of either of these two particle species (and their antiparticles) will lead to much smaller x 2 /dof values while 
the resulting best fit parameters are adjusted within the errors of the parameters resulting from a fit to the full data 
set. In general, the statistical hadronization model tends to systematically over estimate the f2 + f2 rapidity density 
and under estimate the other hyperon yields at 200A GeV. 

Fits to the combined STAR+PHENIX+BRAHMS data shed further light to the issue of short lever arm of the 
PHENIX and BRAHMS data sets. Namely, we have found out that the resulting cross normalization factors fp 
and fs are unrealistically small (« 0.8) if one takes into account the nucleons measured by PHENIX and BRAHMS 
collaborations. Both of these factors can be determined directly from the data by dividing the experimental rapidity 
densities of PHENIX and BRAHMS hadrons by the corresponding ones from STAR collaboration. We have plotted the 
ratios of average pions and kaons, i.e. n + + n~ and K + + K~ from PHENIX divided by the corresponding quantities 
from STAR in fig. \E\ As one can see, the experimental ratios mostly lie between 0.9 and 1.0 at all centralities. 
However, the deviation from unity is clearly large enough so that the previously introduced fp must be implemented. 
In the same figure, we also show the resulting f P (thick line) fitted to the STAR+PHENIX+BRAHMS data when 
excluding protons and anti- protons from PHENIX and BRAHMS. One can see that the fitted fp closely follows the 
experimental ratio of PHENIX and STAR pion rapidity densities. Thus, it seems that the cross-normalization factor for 
the nucleons would be different and much smaller (0.8 or below) than the cross normalization factor for the pions and 
kaons. Unfortunately, this cannot be estimated directly from the data, because, unlike for BRAHMS and PHENIX, 
nucleons from STAR collaboration include all the weak decay products of hyperons. It should be mentioned that both 
BRAHMS and PHENIX weak feeding corrections at 200 A GeV are based on the same PHENIX A and A measurement 
at 130 A GeV. If the underlying assumption of A/p being constant at all centralities at all beam energies at RHIC was 
not correct, then the weak decay corrected proton and anti-proton rapidity densities quoted by PHENIX and BRAHMS 
would be incorrect as well, which could partly explain the failure of statistical hadronization model fits to the PHENIX 
and BRAHMS data alone as well as the unrealistically low fp and fp from fits to the whole data set. In fact, PHENIX 
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(and BRAHMS) weak decay correction is based on the A mc /p exc '=0.89±0.07 [33[ ratio in which A mc denotes the 
inclusive rapidity density of A including weak feeding from S's and f2 while p excl = p — 0.64A mc is the rapidity density 
of protons from which feeding from inclusive A's has been subtracted. This can be compared with the corresponding 
ratio at 200A GeV measured by the STAR collaboration. In order to estimate the inclusive A rapidity density, we sum 
up all the contributions A inc = A excl + (~- -> A) + (0 -> A) + (S° -> A) = 16.7 + 2.17 + 0.68x0.53/2.0 + 0(2). The 3° 
rapidity density is not measured but is expected to be 2° < E~. Thus, we get \ mc /p excl ~ 21/(34.7 — 0.64 x 21) rss 1, 
clearly different from the value obtained by PHENIX collaboration at 130A GeV. For comparison, we note from the 
table U that the STAR data at 130A GeV suggests that again, A inc /p excl w 16.2/(26.4- 0.64 x 16.2) w 1. Because 
of this significant discrepancy in the A/p among the three experiments, we have decided to exclude the proton and 
anti-proton rapidity densities measured by PHENIX and BRAHMS collaborations in the analysis and all results in 
this paper are evaluated excluding these 4 measurements. This way, our fitted cross-normalization factors fp and fg 
follow the actual ratios of pions and kaons determined from the data and provide more reliable estimate compared 
with fits including the p and p from PHENIX and BRAHMS in which cases the low A/p ratio would bias the fit 
towards lower temperatures as well as lower fp and fs- 

We have performed fits to the formula (p} in 7 out of 8 and in 8 out of 11 most central centrality bins at 130 
and 200A GeV respectively. There are two reasons why we have refrained from estimating the freeze-out parameters 
in very peripheral bins. First, the interpolations of lower multiplicity particles, such as hyperons, become less and 
less accurate going to more peripheral collisions and using only light mesons and nucleons makes the fit unstable, 
as has been already discussed. Secondly, in extreme peripheral collisions, the role played by exact conservation laws 
(so-called canonical suppression), especially for multi-strange baryons, may become important. Yet, it is very difficult 
to make a definite assessment of this effect onto rapidity densities rather than full phase space yields. We remind that 
at top SPS energy, the difference in fi's fully integrated yield calculated in grand-canonical and S-canonical (enforcing 
vanishing net strangeness) ensembles is 32% and 14% with Np 16 and 40 respectively |19j. Most likely, these figures do 
not change significantly at RHIC and so it is safer to use the simple grand-canonical formula (|TJ) only when Np > 50. 

IV. COMPARISON WITH PREVIOUS ANALYSES 

Up to date several statistical hadronization analyses similar to ours have been carried out on RHIC Au-Au collisions. 
As discussed in the introduction, the majority of them formed ratios of rapidity densities a posteriori, thus possibly 
introducing a bias in the estimation [2pj ] . Also, some of the data we have been using was not yet available. It is then 
useful to compare our results with previous ones in order to see how much the different input can affect the final 
result. 

In [3, different combinations of ratios of rapidity densities measured by STAR, PHENIX, PHOBOS and BRAHMS 
collaborations in Au-Au collisions at 130A GeV in 3 centrality bins are analyzed. Depending on the used set of ratios 
of rapidity densities, the authors find somewhat different chemical freeze-out conditions. In the case of maximal 
amount of ratios fitted in the analysis, their chemical freeze-out parameters agree with ours within the errors. 

STAR collaboration has determined [H, E3 the statistical hadronization model parameters in Au-Au collisions at 
200 A GeV by using different combinations of ratios of rapidity densities as a function of centrality. Their most recent 
results are in very good agreement with ours. 

The inhomogeneous chemical freeze-out model [l5| takes into account possible fluctuations in the temperature and 
baryon number among created clusters in a collision event. The model has been applied to determine the temperature 
and baryon chemical potential in Au-Au collisions at 130 and 200A GeV at RHIC. In this analysis, ratios of rapidity 
densities are implemented and the authors have found out that the choice of ratios included in the analysis can 
indeed bias the outcome. To try to minimize the bias, particle/antiparticle ratios are included along ratios of various 
particles and negative pions. The central values of the distributions of temperature and baryon chemical potential 
are in approximate agreement with our findings. 

In [lj| ratios of rapidity densities at RHIC energies have been used and the effect of including different sets of 
rapidity densities as well as ratios of them is explored. Admittedly, the authors find a difference in the fit outcome 
whether using rapidity densities or ratios or different set of ratios. Among the many quoted results, some are in 
agreement with ours. 

To our knowledge, the only other statistical model analysis [13j which used rapidity densities themselves instead 
of forming ratios includes the PHENIX ir^ , p and p as well as two ratios measured by the STAR collaboration, 
and finds somewhat lower temperatures than we do, for the same version of the statistical model. We deem that this 
discrepancy owes to the poorer data set available when that analysis was carried out. 

In conclusion, several groups have analyzed the rapidity densities at RHIC 130 and 200 A GeV and the results of 
all groups seem to agree rather well. The resulting baryon chemical potentials agree very well among all groups and 
seem to be fairly insensitive both to the set of particle species included in the analysis as well as to the details of 
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the version of the statistical hadronization model. Chemical freeze-out temperature, on the other hand, shows larger 
discrepancies and seems to be more sensitive to the input data set, which, as has been emphasized, is an effect to be 
expected in fitting different subsets of ratios, also without including correlations. Finally, the values and behaviour 
of 75 as a function of centrality are in very good agreement with previous findings, especially with the ones fl4l . Il7j 
calculated with the THERMUS [H[ package. 

On the other hand, we observe a worse fit quality at 200A GeV than generally reported by previous analyses. This 
difference is certainly owing to our updated and more accurate data set but, again, to some extent this is possibly 
related to having fitted rapidity densities instead of a subset of ratios. 




50 100 150 200 250 300 350 400 



FIG. 5: The average rapidity density of pions (open circles) and kaons (open squares) measured by PHENIX collaboration 
divided by the same quantity measured by STAR in Au-Au collisions at 200A GeV as a function of centrality. The fitted 
cross-normalization factor fp (thick line) closely follows the measured ratio of pions. Dashed lines visualize the errorband of 
fp- 



V. SUMMARY AND CONCLUSIONS 



In summary, we have analyzed within the statistical hadronization model the rapidity densities of various hadron 
species at mid-rapidity in Au-Au collisions at 130 and 200A GeV measured by STAR, PHENIX and BRAHMS 
collaborations and determined the relevant statistical hadronization model best fit parameters. This completes our 
previous analyses at lower center-of-mass energies measured at SPS and AGS. 

We have used as input data for the analysis only rapidity densities and not ratios formed out of them because of 
the bias introduced in fitting subset of ratios. Although a direct comparison is not possible because the data set used 
in this analysis is the most up-to-date, our results are in general good agreement with those of previous analyses, 
showing that the effective value of the bias introduced by the actual choice of ratios therein was small and most likely 
within the fit error. 

We have found out that the statistical hadronization model, as implemented by the formula |T]) can describe the 
rapidity densities measured at RHIC relatively well, although discrepancies between data and model are visible and 
larger than some other groups using ratios in their analyses have reported. This is also reflected in the higher x 2 /dof 
value that we find at V*nn = 200 GeV compared to that at \/i NN = -^0 GeV. Since the relative deviations between 
data and model are comparable at both energies, we conclude that the better relative accuracy of measurements at 
the higher energy has overcome the theoretical accuracy of the simple formula (p}. We interpret this not as a failure 
of the statistical model itself, but an indication that corrections to the simple assumptions underlying formula ([T]) 
would be necessary, like those discussed in Sect. IIII1 
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FIG. 6: Chemical freeze-out temperature as a function of center-of-mass energy per participant pair (left panel (a)) and as a 
function of baryon chemical potential (right panel (b)) in central heavy-ion collisions at AGS, SPS and RHIC. The chemical 
freeze-out shown in the right panel as well as the curve shown in the left panel are empirical fits to the AGS and SPS points 
taken from our previous paper [l9| ]. 



A major result of our analysis is the stability of the temperature as a function of centrality, especially at v^nn = ^00 
GeV, where all values range from about 166 to 171.4 MeV, hence with an overall spread of around 3%. This confirms 
previous findings from STAR collaboration. 

The strangeness under-saturation parameter increases mildly from peripheral to central collisions where it almost 
attains 1. Therefore, RHIC data in peripheral collisions demonstrates the phenomenon of phase space under-saturation 
for mid-rapidity yields. This is also observed in the dependence of normalized <f> meson yield as a function of centrality 
|36j |: since the temperature is essentially constant, this behaviour can only be parameterized with a 75 varying as 
a function of centrality. The authors have recently proposed [44j an explanation of strangeness under-saturation in 
terms of a superposition of NN collisions and a completely equilibrated hadronic system originated from the central 
core of the collision, where the Quark Gluon Plasma is formed. Such a scenario will be investigated in more detail in 
future works. 

Finally, we find that the fitted chemical freeze-out temperatures and baryon chemical potentials in central Au-Au 
collisions nicely fit previously extrapolated curves from lower heavy ion collision energies as shown in fig. [SJ 
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TABLE I: Our estimates for the STAR (top panel) and PHENIX (middle panel) rapidity densities of various hadron species 
in different centrality windows in Au-Au collisions at 130^4 GeV. STAR K ± , p and p are experimental values while all others 
are derived from the measured values shown bold face in the bottom part of the table by interpolation described in Section 
Iff Al The STAR pion rapidity densities are corrected for weak decays while all other rapidity densities include the weak decay 
products (if any) of weakly decaying resonances. 
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centrality 


0-5% 


5-10% 


10-15% 


10-20% 


15-20% 




N P 


351.4±2.9 


299.0±3.8 


253.9±4.3 


235±9 


215.3±5.3 




n+ [12, 45] 


322.2±19.2 


257.0±15.2 


210.8±12.7 


193.8±11.4 


176.6±10.7 




7T- [12, 45] 


327.0±19.5 


260.7±15.4 


213.7±12.8 


196.1±11.6 


178.7±10.8 




K+ [12, 45] 


51.27±5.92 


40.82±4.25 


32.9±3.2 


29.97±2.86 


27.0±2.6 




K- [12, 45] 


49.47±5.71 


39.78±4.15 


31.7±3.1 


28.74±2.74 


25.8±2.4 


s 


P [12, 45] 


34.70±4.10 


28.23±2.99 


22.0±2.2 


20.12±1.94 


18.3±1.7 


T 


p [12, 45] 


26.70±3.15 


21.42±2.27 


17.1±1.7 


15.69±1.51 


14.3±1.4 


A 


A [17] 


16.7±1.1 


13.55±0.91 


11.02±0.77 


10.0 ±0.7 


8.98±0.64 


R 


A [17] 


12.7±0.9 


10.36±0.69 


8.47±0.56 


7.7 ±0.5 


6.93±0.46 




B- [lJJ 


2.17±0.20 


1.86±0.14 


1.55±0.11 


1.41±0.09 


1.268±0.074 




E+ [17] 


1.83±0.21 


1.52±0.14 


1.25±0.11 


1.14±0.09 


1.027±0.075 




</> [35 , 36 , 45J 


7.95±0.74 


6.81±0.73 


5.89±0.53 


5.37±0.51 


4.82±0.51 




n + H [17] 


0.53±0.06 


0.445±0.044 


0.368±0.035 




0.299±0.030 


P 


7T+ [34] 


286.4±24.2 


239.6±20.5 


204.6±18.0 




173.8±15.6 


H 


7T" [34] 


281.8±22.8 


238.9±19.8 


198.2±16.7 




167.4±14.4 


E 


if+ [34] 


48.9±6.3 


40.1±5.1 


33.7±4.3 




27.9±3.6 


N 


if~ [34] 


45.7±5.2 


37.8±4.3 


31.1±3.5 




25.8±2.9 


I 


P [34] 


18.4±2.6 


15.3±2.1 


12.8±1.8 




10.6±1.5 


X 


P [34] 


13.5±1.8 


11.4±1.5 


9.5±1.3 




7.9±1.1 


B 


7T+ 


309.8±32.3 


256.8±25.2 


213.7±21.3 


196.2±19.7 


178.8±18.2 


R 


7T~ 


302.6±31.6 


253.2±24.9 


212.2±21.2 


195.3±19.6 


178.4±18.2 


A 


K+ 


49.8±5.2 


40.2±4.0 


32.7±3.3 


29.8±3.0 


26.9±2.8 


H 


K~ 


44.7±4.7 


37.1±3.7 


30.7±3.1 


28.1±2.9 


25.5±2.6 


M 


P 


20.0±2.1 


17.2±1.7 


14.6±1.5 


13.4±1.4 


12.2±1.3 


S 


P 


14.8±1.6 


12.6±1.2 


10.6±1.1 


9.7±1.0 


8.8±0.9 




centrality 


0-10% 










N P 


328±6 








B 


7T+ [3JJ 


283.3±28.4 








R 


7T" [37] 


277.9±27.9 








A 


if+ [37] 


45.0±4.55 








H 


if" [37] 


40.9±4.14 








M 


P [37] 


18.6±1.87 








S 


P [37] 


13.7±1.38 









TABLE II: Rapidity densities of various hadrons in Au-Au collisions at 200A GeV in different centrality windows. Bold 
face numbers denote measured values while numbers written with standard fonts denote our estimates. Our estimates are 
interpolated from the experimental values as described in the Sect III Al The STAR p and p rapidity densities do include weak 
feeding from (multi)strange hyperons while PHENIX and BRAHMS p and p rapidity densities include feeding from E's only. 
All other rapidity densities in this table are corrected for the weak feeding (if relevant). 
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centrality 


20-30% 


30-40% 


40-50% 


50-60% 


60-70% 


70-80% 




N P 


166.6±5.4 


114.2±4.4 


74.4±3.8 


45.5±3.3 


25.7±3.8 


13.4±3.0 




7T+ [12, 45] 


134.93±7.78 


89.24±5.13 


58.66±3.35 


36.24±2.07 


21.07±1.20 






7T- [12, 45] 


136.07±7.84 


89.64±5.16 


58.85±3.36 


36.33±2.07 


21.13±1.20 






K+ [12, 45] 


20.48±1.77 


13.61±1.11 


8.690±0.680 


5.400±0.410 


2.980±0.220 






K~ [12, 45] 


19.68±1.70 


13.18±1.07 


8.370±0.660 


5.190±0.390 


2.890±0.220 




s 


P [12, 45] 


14.39±1.26 


9.300±0.760 


6.170±0.480 


3.880±0.290 


2.200±0.160 




T 


P [12, 45] 


11.180±0.980 


7.460±0.610 


4.930±0.390 


3.160±0.240 


1.840±0.140 




A 


A 


6.67±0.47 


4.39±0.34 


271±0.21 


1.42±0.16 






R 


A 


5.18±0.35 


3.43±0.27 


2.13±0.17 


1.15±0.12 








H — 


0.903±0.037 


0.537±0.031 


0.315±0.025 


0.205±0.035 








H 


0756±0.043 


0.484±0.037 


0.295±0.029 


0.165±0.031 








<t> [35, 36, 45] 


3.47±0.44 


2.29±0.23 


1.44±0.14 


0.810±0.092 


0.450±0.051 


0.20±0.022 




n + n 


0.213±0.025 


0.127±0.020 


0.073±0.011 


0.053±0.015 






P 


7T + [34] 


130.3±12.4 


87.0±8.6 


54.9±5.6 


32.4±3.4 


17.0±1.8 


7.9±0.8 


H 


7T- [34] 


127.3±11.6 


84.4±8.0 


52.9±5.2 


31.3±3.1 


16.3±1.6 


7.7±0.7 


E 


[34] 


20.6±2.6 


13.2±1.7 


8.0±0.8 


4.5±0.4 


2.2±0.2 


0.89±0.09 


N 


if - [34] 


19.1±2.2 


12.3±1.4 


7.4±0.6 


4.1±0.4 


2.0±0.1 


0.88±0.09 


I 


P [34J 


8.1±1.1 


5.3±0.7 


3.2±0.5 


1.8±0.3 


0.93±0.15 


0.40±0.07 


X 


P[34J 


5.9±0.8 


3.9±0.5 


2.4±0.3 


1.4±0.2 


0.71±0.12 


0.29±0.05 


B 


7T+ 


139.2±13.9 


100.2±11.4 


67.1±9.3 


27.9±77 






R 


7T~ 


138.9±13.9 


98.7±11.2 


64.0±9.1 


27.7±7.7 






A 


AT+ 


20.8±2.1 


15.1±1.8 


10.0±1.5 


2.7±1.3 






H 


if- 


19.4±2.0 


13.3±1.6 


8.4±1.4 


3.4±1.2 






M 


P 


9.22±0.93 


5.98±071 


3.61±0.58 


2.02±0.54 






S 


P 


6.70±0.68 


4.50±0.53 


2.86±0.46 


1.58±0.42 








centrality 


20-40% 


40-60% 


60-80% 




N P 


141±8 


62±9 


21±6 


S 


A [17] 


5.53±0.39 


2.07±0.14 


0.58±0.04 


T 


A [17] 


4.30±0.30 


1.64±0.11 


0.48±0.03 


A 


S" [17] 


0.72±0.03 


0.26±0.02 


0.063±0.005 


R 


E+ [IT] 


0.62±0.04 


0.23±0.02 


0.061±0.004 




Q + Q [171 


0.17±0.02 


0.063±0.009 






B 


7T+ [37] 


119.7±12.1 


47.5±4.85 






R 


7T- [37] 


118.8±12.0 


46.3±4.71 






A 


*f+ [37] 


17.9±1.83 


6.3±0.68 






H 


if- [37] 


16.3±1.7 


5.9±0.64 






M 


P [37] 


7.6±0.77 


2.81±0.29 






S 


P [37] 


5.6±0.6 


2.22±0.24 







TABLE III: Rapidity densities of various hadrons in Au-Au collisions at 200A GeV in different centrality windows. Bold 
face numbers denote measured values while numbers written with standard fonts denote our estimates. Our estimates are 
interpolated from the experimental values as described in the Sect III Al The weak decay corrections are the same as listed in 
the previous table. 
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centrality 


T [MeV] 


Ms [MeV] 


7s 


dV/dy [fm 3 ] 


X 2 /dof 


fp 




STAR 130A GeV best fit parameters 


0-6% 


165.9±5.3 


35.1±12.6 


1.109±0.078 


1097±258 


5.9 / 7 




6-11% 


165.3±5.1 


34.0±12.2 


1.104±0.075 


925±211 


5.4 / 7 




11-18% 


165.3±5.3 


33.2±12.6 


1.053±0.075 


760±176 


6.1 / 7 




18-26% 


162.5±5.4 


30.7±13.3 


0.977±0.078 


712±173 


2.6 / 7 




26-34% 


163.1±2.5 


28.1±15.0 


0.907±0.056 


534±61 


1.4 / 7 




34-45% 


161.1±7.4 


27.3±21.8 


0.863±0.091 


410±140 


2.5 / 7 




45-58% 


153.5±7.8 


26.9±26.1 


0.823±0.096 


352±136 


4.7 / 7 






PHENIX 130A GeV best fit parameters 


0-5% 


161.4±6.9 


33.9±16.8 


1.02±0.16 


1397±475 


1.7 / 2 






STAR+PHENIX 130A GeV best fit parameters 


0-6% 


163.8 ±4.1 


36.3 ±10.2 


1.109 ±0.067 


1225 ±228 


8.2 / 12 


0.919±0.067 


6-11% 


163.7 ±4.0 


36.1 ±9.91 


1.097 ±0.064 


1013 ±182 


8.6 / 12 


0.911±0.065 


11-18% 


163.8 ±4.0 


35.9 ±10.1 


1.043 ±0.064 


833 ±150 


10.3 / 12 


0.897±0.065 


18-26% 


161.9 ±3.9 


34.4 ±10.7 


0.954 ±0.024 


750 ±133 


7.2 / 12 


0.844±0.060 


26-34% 


162.1 ±4.4 


32.9 ±5.26 


0.884 ±0.067 


576 ±116 


7.1 / 12 


0.801±0.061 


34-45% 


159.1 ±4.9 


32.7 ±13.6 


0.827 ±0.071 


460 ±106 


7.5 / 12 


0.798±0.064 


45-58% 


153.3 ±5.4 


26.8 ±18.3 


0.813 ±0.080 


357 ±96 


5.0 / 12 


0.841±0.018 



TABLE IV: Statistical hadronization model best fit parameters at chemical freeze-out in Au-Au collisions at 130A GeV. 



centrality 


T [MeV] 


/is [MeV] 


7s 


dV/dy [fm 3 ] 


X 2 /dof 


fp 


fB 




STAR 200,4 GeV best fit parameters 


0-5% 


168.0±6.2 


28.8±14.3 


0.935±0.064 


1419±377 


22.2 / 8 






5-10% 


169.5±6.8 


29.1±14.7 


0.941±0.069 


1055±304 


26.8 / 8 






10-15% 


167.0±7.3 


26.8±16.2 


0.979±0.079 


941±296 


36.9 / 8 






10-20% 


168.8±6.0 


27.0±13.3 


1.054±0.077 


745±190 


22.4 / 7 






15-20% 


167.8±7.0 


27.1±15.3 


0.971±0.076 


750±225 


36.1 / 8 






20-30% 


169.2±6.3 


27.9±12.5 


0.954±0.064 


537±141 


30.0 / 8 






30-40% 


166.4±5.5 


22.0±12.7 


0.951±0.063 


399±95 


26.7 / 8 






40-50% 


165.8±5.2 


21.4±12.8 


0.900±0.059 


274±62 


25.1 / 8 






50-60% 


164.9±3.1 


20.5±8.0 


0.902±0.043 


173±24 


7.0 / 8 








STAR±PHENIX±BRAHMS 200A GeV best fit parameters 


0-5% 


169.2±5.2 


29.5±11.2 


0.929±0.044 


1336±302 


23.4 / 14 


0.863±0.073 


0.89±0.12 


5-10% 


171.2±5.2 


29.7±11.3 


0.928±0.048 


976±210 


28.1 / 14 


0.912±0.078 


0.931±0.082 


10-15% 


168.9±5.5 


27.6±12.6 


0.960±0.054 


868±201 


39.7 / 14 


0.928±0.093 


0.932±0.096 


15-20% 


169.9±5.4 


27.9±12.0 


0.951±0.051 


686±156 


38.9 / 14 


0.948±0.095 


0.944±0.098 


20-30% 


171.4±5.0 


28.8±9.90 


0.935±0.046 


489±102 


32.8 / 14 


0.942±0.088 


0.970±0.091 


30-40% 


168.3±4.4 


23.1±10.2 


0.930±0.043 


367±68 


30.0 / 14 


0.931±0.085 


1.04±0.10 


40-50% 


167.3±3.1 


22.5±9.35 


0.883±0.030 


257±34 


27.9 / 14 


0.870±0.068 


1.06±0.12 


50-60% 


166.2±3.3 


20.8±10.4 


0.874±0.045 


165±24 


11.1 / 14 


0.818±0.051 


0.67±0.11 



TABLE V: Statistical hadronization model best fit parameters at chemical freeze-out in Au-Au collisions at 200A GeV. Errors 
are scaled according to the PDG scaling scheme [411 ] by a factor ^x 2 /dof. 
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Experiment (E) Model (M) Residual (M - E)/E (%) 



STAR Au-Au 130A GeV 0-6% most central collisions 



7T + 


234±24 


213 


-0.87 


-8.93 




234±24 


215 


-0.76 


-7.82 


K+ 


46.2±6.0 


49.0 


0.46 


5.96 




41.9±5.4 


45.7 


0.69 


8.95 


V 


26.4±5.8 


32.2 


1.0 


22.3 


V 


18.7±4.1 


22.4 


0.90 


19.8 


4> 


6.26±0.90 


6.61 


0.38 


5.48 


A 


16.2±1.6 


16.3 


0.064 


0.648 


A 


11.8±1.2 


12.2 


0.32 


3.28 




2.18±0.28 


2.00 


-0.63 


-8.24 




1.87±0.24 


1.59 


-1.2 


-14.9 



TABLE VI: Comparison between estimated rapidity densities in the combined fit and rapidity densities measured by STAR 
in central Au-Au collisions at 130A GeV. The 3rd and 4th column show the discrepancy between data and model in units of 
standard error and in percentages, respectively. 



Experiment (E) Model (M) Residual (M - E) /E (%) 



STAR Au-Au 200A GeV 0-5% most central collisions 



7T + 


322±19 


326 


0.21 


1.24 




327±20 


328 


0.055 


0.326 


K+ 


51.3±5.9 


57.1 


0.99 


11.4 


K~ 


49.5±5.7 


53.9 


0.78 


8.97 


V 


34.7±4.1 


41.8 


1.7 


20.4 


V 


26.7±3.1 


30.9 


1.3 


15.9 


<t> 


7.95±0.74 


6.73 


-1.6 


-15.3 


A 


16.7±1.1 


14.4 


-2.1 


-13.9 


A 


12.70±0.92 


11.07 


-1.8 


-12.8 




2.17±0.20 


2.02 


-0.74 


-6.86 




1.83±0.21 


1.67 


-0.76 


-8.51 


n + n 


0.530±0.057 


0.651 


2.1 


22.9 



TABLE VII: Comparison between estimated rapidity densities in the combined fit and rapidity densities measured by STAR 
in central Au-Au collisions at 200A GeV. The 3rd and 4th column show the discrepancy between data and model in units of 
standard error and in percentages, respectively. 



